Effects of electron coupling to intra- and inter-molecular vibrational modes on the 
transport properties of single crystal organic semiconductors 
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Electron coupling to intra- and inter-molecular vibrational modes is investigated in models ap- 
propriate to single crystal organic semiconductors, such as oligoacenes. Focus is on spectral and 
transport properties of these systems beyond perturbative approaches. The interplay between dif- 
ferent couplings strongly affects the temperature band renormalization that is the result of a subtle 
equilibrium between opposite tendencies: band narrowing due to interaction with local modes, band 
widening due to electron coupling to non local modes. The model provides an accurate description 
of the mobility as function of temperature: indeed, it has the correct order of magnitude, at low 
temperatures, it scales as a power-law T~ s with the exponent 8 larger than unity, and, at high 
temperatures, shows an hopping behavior with a small activation energy. 



I. INTRODUCTION 

In the last years, the availability of single crystal 
organic field-effect (OFET) transistors has represented 
a step forward in the field of plastic electronics since 
these systems show a charge mobility at least one or- 
der of magnitude larger than that of thin filmsj- Among 
them, the most promising and studied are those based 
on oligoacenes, such as pentacene and rubrene3 

Although the transport properties in the linear regime 
have been carefully extracted by studying clean single 
crystal OFET devices, the intrinsic transport mechanism 
acting in organic semiconductors (OS) is, up to now, not 
fully understood. Typically, in these materials, the in- 
duced doping is not very high (typically much smaller 
than one charge carrier for ten molecules), so that the 
mobility is analyzed in transport measurements. At tem- 
peratures close or higher than 100-ftf, the mobility /i p of 
these systems exhibits a power-law behavior ~ T , 
with S ~ 2), 2-4 which cannot be simply ascribed to band 
. transport £ Moreover, in some systems, starting from 



room temperature, there is a crossover from band-like 
to activated hopping behavior.— For example, in naph- 
thalene and anthracene, along the a and b axis, there is 
a change in the temperature behavior of mobility, and, 
along the c axis, an activated regime takes place with a 
small activated gap of the order of 20meV£ Moreover, in 
pentacene, above room temperature, the mobility shows 
an upturn with increasing temperature^ 

Recently, many photoemission experiments have been 
performed in pentacene and rubrene.— ~— The quasi- 
particle energy dispersion does not exhibit a strong mass 
renormalization indicating moderate values of electron- 
phonon (el-ph) coupling. For pentacene, the band struc- 
ture is reduced only by about 15% going from 75K to 
300K. 

Typically, the crossover from tunneling to hopping be- 
havior observed in transport properties is ascribed to 
the formation of local polaronS ) 12 i 13 and described with 
the simplest model for molecular crystals, the Holstein 
modeh 14 ' 15 The modeling assumes that the most rele- 
vant carrier interaction is with local modes whose fre- 



quency is lower than the typical hopping energy i 12 ' 13 
The used approaches only generalize the previous treat- 
ment of Holstein works / 4 ' 15 starting from the question- 
able limit of narrow electronic band even for the coupling 
to low-frequency modes. Actually, recent results, based 
on the adiabatic approach for the vibrational modes, have 
pointed out that, up to room temperature, within a pa- 
rameter range valid for many oligoacenes, polaronic ef- 
fects do not play a prominent role if low frequency modes 
are excited in the Holstein modeli 1 ^ 

Recent ab-initio calculations have evidenced that 
charge carriers are mostly coupled to intra-molecular 
modes with high frequency. Moreover, it has been 
found that another very relevant interaction is with inter- 
molecular modes with low frequency in comparison with 
typical electron hoppings£ii£ Another important result 
is that the reorganization energy (related to the polaron 
binding energy) decreases with increasing the number of 
benzene rings in oligoacenes (for example, going from 
naphthalene to pentacene). 

In view of ab-initio results, Holstein-Peierls models 
have been considered for the description of OSjiSr— mod- 
elling the carrier interaction with low-frequency inter- 
molecular modes via a modulation of the transfer hop- 
ping. Three dimensional models have been mainly 
applied to naphtalenei 1 ^ One-dimensional models sim- 
ilar to Su-Schrieffer-Heeger (SSH) model™ have been 
introduced to study the OS along the relevant stack 
direction^ 2 - - — These models treat the inter-molecular 
modes as classical exploiting their low energy in compar- 
ison with the electronic hopping scale. The approaches 
are non perturbative in the electron-phonon coupling. 
The mobility is quite well described at low temperatures 
showing a robust power-law T~ s . However, the power- 
law behavior extends up to high T and there is no sign 
of upturn in the mobility as function of temperature. 

Aim of this work is to clarify the effects of electron 
coupling to intra- and inter-molecular vibrational modes 
on the transport properties of prototype single crystal 
OS, such as oligoacenes. To this aim, a one-dimensional 
model with coupling to intra- and inter-molecular modes 
has been analyzed. We will show that the interplay be- 
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tween local and non local el-ph interactions is able to 
provide a very accurate description of the mobility and 
to shed light on the intricate mechanism of band narrow- 
ing with increasing temperature. 



Finally, in Eq. ([T]), H e i-i nter represents the SSH-like 
term with electron coupling to inter-molecular modes 

Hel-Inter = O/nter ^^iVi+l—Vi) fcfcj+1 + 4+l C i) > ( 5 ) 



II. MODEL 

We consider a one-dimensional model with coupling to 
intra- and inter-molecular modes similar to one recently 
introduced, where the treatment only concerns the study 
of spectral properties." 2 - 6 - The coupling to intra-molecular 
modes is Holstein-like, that to inter-molecular modes is 
SSH-like. It can be summarized in the following hamil- 
tonian: 

H = H el ) + H Intra + H Inter + H el-Intra + H e i- Inter . (1) 

(0) 



In Eq. (J]) , the free electronic part is 



H 



(0) 



c !+i c -' 



(2) 



where t is the bare electron hopping between the near- 
est neighbors on the chain, c\ and Ci are the charge 
carrier creation and annihilation operators, respectively, 
relative to the site i of a chain with lattice parame- 
ter a. For the transfer hopping the ab-initio estimate 
is: t ~ 50 — 100meV\— We consider a single-band one- 
dimensional electronic structure since it represents the 
simplest effective model in the anisotropic OS to ana- 
lyze the low energy features responsible for the mobility 
properties. 

In Eq. ([T]), Ha \ with a = Intra, Inter is the Hamil- 
tonian of the free optical molecular modes 



k 7 2 



(3) 



where z a ,i and r a ^i are the oscillator displacement and 
momentum of the mode a, respectively, m a the oscillator 
mass and k a the elastic constant of the mode a. For non- 
local modes, we fix = zi nter .i, for local modes, Xi = 
zintra,i and Pi — rj ntra ,i- The inter- molecular modes are 
characterized by small frequencies (tvjji nte r — 5 — WmeV) 
in comparison with the transfer hopping. 6 ' 22 On the con- 
trary, the most coupled intra-molecular modes have large 
frequencies (fiw/nt™ — 130 — \80meV)& 

In Eq. (J]), H e i~ intra is the Holstein-like Hamiltonian 
describing the electron coupling to intra-molecular modes 



H e i- 
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Otlntr 
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with ai n tra coupling constant to local modes and ni — 
c\ci local density operator. The dimensionless constant 
9 intra = ct Intr a/ y / 2hm Intra uj^ ntra is used to describe this 
el-ph coupling,— In single crystal OS, g intra is in the 
weak to intermediate regime (of the order of unity). 



with ctinter coupling constant to non local modes. In 
the adiabatic regime for non local modes (fiw/„ ter <C t), 
the dimensionless quantity Xinter — ctj nter /4ki n t e rt fully 
provides the strength of the electron coupling to inter- 
molecular modes. The typical values of A are in the in- 
termediate coupling regime (of the order of or less than 

o.i). 

In the following, we will use units such that lattice pa- 
rameter a = 1, Planck constant K = 1, Boltzmann con- 
stant ks = 1) and electron charge e = 1. We will analyze 
systems with chain size L = 32 (where the system reaches 
the thermodynamic limit) and we will measure energies 
in units of t ~ 80meV. We fix lj intra — 2.0 as the highest 
energy scaled For the calculation of dynamic quantities, 
an additional small broadening T is introduced in order 
to simulate the effect of a tiny disorder. We fix T — 0.05<, 
therefore it is the smallest energy scale. The results at 
low temperatures are only slightly dependent on the value 

of r. 



III. METHOD 

Since a very low carrier density is injected into OS, we 
will study the case of a single particle. 

The temperature range where intrinsic effects are rele- 
vant is w inter < T <C < < u> intra- Therefore, the dynam- 
ics of intermolecular modes is assumed classical. Actu- 
ally, the electron dynamics is strongly influenced by the 
statistical "off diagonal" disorder, that, in the limit of low 
carrier density, is described by the probability function 
P ({yi}) of free classical harmonic oscillators: 



P({Vi}) = 



f/3k 



Inter 



V 2tt 



L/2 



exp 



n Winter \ ^ 9 



. (6) 



with /3 = 1/T. 

At a fixed configuration of non local displacements 
{yi}, Eq. (JlJ is equivalent to a Holstein model, where the 
electron hopping between nearest neighbor sites is not ho- 
mogeneous and depend on the specific pair (see Eq.©). 
The resulting inhomogeneous Holstein model can be ac- 
curately studied within the modified Lang-Firsov ap- 
proach via a unitary transformation U, provided that the 
system is in the anti-adiabatic regime {u intra > t), where 
the quantum nature of phonons cannot be neglectedj 27 i 28 
The electron is renormalized by the coupling with local 
modes increasing its mass, and the oscillators are dis- 
placed from their equilibrium position to a distance pro- 
portional to the el-ph interaction. In our case, the new 
Hamiltonian H = XJ^HXJ is built through the varia- 



3 



tional transformation: 



U ({yj}) = exp 



^9lntra 



Wllntra^Intra 



E^ ({yj})pi 



where the variational parameters fi ({yj}), giving the 
new centers for the local oscillators at the sites i, have to 
be determined for each fixed configuration of non local 
displacements {yj}. Within our approach ; 28 ' 29 the solu- 
tion of the full problem can be obtained by minimizing 
the free energy of the following inhomogeneous effective 
hamiltonian: 



H eff ({ Vl }) = tfjl a 



E 



[fi ({%» - 2/, ({%})] rn 



i 

where 



-t + OtlnteriVi+l ~ Vi)] e T '-'+ 1 «»» (c\c l+1 + /i.e.) , 



Tij ({Vj}) = 



9 2 mtra [ff ({%}) + //({%})] (2^Vo + l) 



(8) 

represents the term giving the reduction of the bond hop- 
ping, with No = 1/ '( y e^ Ulntra — 1) Bosc distribution of the 
local modes. The distribution of fi values depends on 
the distribution P ({yi}). For example, close to T = 0.5t, 
the deviations of /, from the average value can become 
of the order of twenty per cent in the parameter regime 
investigated in this work. For strong el-ph coupling and 
high temperatures, the values of fi are close to one (X^i+i 
large) , meaning that the system is characterized by small 
polaron behavior and an incoherent hopping dynamics. 
For intermediate el-ph coupling and low temperatures, 
the distribution of fi can have an average smaller than 
unity. This implies that there is a possibility of tunnel- 
ing dynamics with a renormalized bandwidth (T^i+i not 
large) . 

At fixed configuration {yj}, one has to diagonalize 
H e ff yielding L eigenvalues E n and eigenvector compo- 
nents bj(n), with j = 1, .., L. In this paper, we will focus 
on the density of states D(uj) ({yj})- 



%)(fe}) = 7E4iM ({yj}), 



(9) 



where A i:i (tu) ({yj}) represents the diagonal term of the 
spectral function 

oo 



with Vi ({%}) = 2gj ntra ff({y J })^/N (N + l), k(z) 
modified Bessel function of order I and the function B n (z) 



1 F 

Bn(z) = ~- 



ir(z-E n ) 2 + rf 



(10) 



The sum over I provides the contribution of phonon repli- 
cas to the density of states.— 

Another central quantity of this work is the mobility 
as function of the temperature. At fixed configuration 
{yj}, the mobility fJ-({yj}) of the inhomogeneous Hol- 
stein model is determined starting from the real part of 
the conductivity Re[a(oj)] ({yj}), taking the limit of zero 
frequency and dividing for the single particle density 1/L 
of the system. In the linear regime, the real part of the 
conductivity is derived from the Kubo formula 2 ^ 



Re[a(w)] ({ yj }) 
where 

3 = 



(1 - e-g") 1 
2lu L 



dte^UHMO)), 



^E^J+^A^ ({yj}) 

j,s 



(11) 

(12) 



is the current operator, hj t s ({yj}) — t~ a ^(yj+s—yj) de- 
notes the generalized hopping and 5 = —1,1 indicates the 
nearest neighbors. In our approach, the current-current 
correlation function becomes 

(i f (i)i(0)) = E T,^ h ^n,s 1 C 3 , Sm j 1 (t)E ]iS . Ju5l (t), 

j,ji S,Si 

where Cj t s-j lt Si (t) is the electron correlation function 
Cj,S; juSl (t) = (c}(t)c J+s (t)c] i+Si (0)c n (0))^ ff , (14) 

and Ej t s-j 1 .8 1 (t) is the phonon correlation function 
Ej,5;jusM = (^(t)X j+s (t)xj i+Si (0)^(0))^, 

with Xi the multi-phonon operator 

fi({yj})p 



Xi = exp 



^Qlntra 



ITlIntra^ Intra 

The function Ej t s-j 1 ,s 1 (t) can be written as 

E j,5;ji,&i(t) 



(15) 
(16) 



e T 3,i+«({y 3 }) e T ii,ii+«i({Vj}) _|_ 



With the first term of the previous equation, one can 
build up the coherent part of conductivity, where the 
charge transfer is not accompanied by processes changing 
the number of phonons. With the second term, retaining 
only the main autocorrelation contribution j = ji, one 
gets the incoherent part of the conductivity, which takes 
into account inelastic scattering processes of emission and 
absorption of phonons ^ 

If an observable O depends on the displacements {yi}, 
first one makes the average over the eigenstates and 
eigenvectors of H e ff({yi}), then over the distribution 
P ({yi}) making the integral 



(0) = 




P({yi})0({ yi }) (17) 
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by means of a Monte-Carlo procedure. Actually, we gen- 
erate a sequence of random numbers distributed accord- 
ing to P {{yi}). For the systems investigated in this pa- 
per, a few thousands of iterations are sufficient to get a 
good accuracy even for dynamic quantities. 

The method exposed above is very accurate in the 
regime uj inter *C t and uj intra > t appropriate to OS. It 
properly takes into account the quantum effects of high 
frequency local vibrational modes. Moreover, the ap- 
proach is able to include spatial correlations relevant in 
quasi one-dimensional systems, in particular vertex cor- 
rections in the calculation of mobility. 



IV. RESULTS 




FIG. 1. Mobility y, as function of temperature for different 
values of local coupling gi n tra- In the inset, the density of 
states D as a function of the frequency at T — for different 
values of local coupling gintra- 

First, we analyze the case Xinter = 0. In Fig. 1, we 
report the mobility with varying the values of g intra- 
For weak gintra, a coherent behavior is present with a 
power-law of the order of 1/T at low temperature. As ex- 
pected, with increasing gintra, there is a crossover from 
tunneling to hopping behavior due to the formation of 
localized small polarons. That occurs at temperatures 
slightly larger than 0.5t for gintra = 13. For large values 
of gintra, the low temperature power-law part is strongly 
reduced and the activated behavior is followed by the 
residual scattering regime at very high 

In the variational approach, the crossover from band- 
like to hopping regime is marked by the value / = I on 
all the sites. However, in the intermediate regime (for 
example gintra — 1-3), at temperatures characterized by 
a metallic behavior, the calculation provides values of / 
smaller than unity (/of the order of 0.6) suggesting that 
the tunneling regime is reduced but still present. More- 



over, at higher temperatures, where / = I is recovered, 
the activated gap A in the hopping mobility is not very 
large: A ~ 0.54*. 

In the inset of Fig. 1, we report the density of states at 
T = for the same values of gintra of the main figure. As 
expected, with increasing gintra, there is a reduction of 
the fundamental polaronic band with a transfer of spec- 
tral weight to high satellite bands appearing at multiples 
of the vibrational frequency u!i n t ra ^ 



m , =2.0, g, . =1.3 

Intra ' °lntra 



T=0.00 
T=0.20 
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X, , =0.09 
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FIG. 2. The density of states as function of the frequency for 
different temperatures at Xinter = 0.09 and gintra = 1.3. 

The description of mobility based only on the intra- 
molecular mode electron coupling fails for two reasons in 
the comparison with low temperature experiments. First, 
it does not provide large values of the mobility (expressed 
in our natural units, that is in terms of /in = ea 2 /h ~ 
7cm 2 /(V-s), taking a = 7l 6 ). For example, at T = 0.314 
(room temperature), fi is less than 0.5/xq for gintra = 13. 
Then, the low temperature mobility due local coupling 
scales as 1/T. 

The next step is to combine the effects of high fre- 
quency local vibrational modes with non local low fre- 
quency ones. These non local modes are considered as 
relevant for the description of the mobility up to room 
temperature. Actually, for rubrene, the carrier mobil- 
ity is dominated by inter-molecular phonons (the non 
local coupling is estimated to be in the intermediate 
regime Xinter — 0.09) since the interaction with the intra- 
molecular modes is almost negligible (gintra — 0.7) M- In 
pentacene, the coupling with local modes is of the same 
order of that in rubrene^ For these materials, our ap- 
proach, being the local el-ph coupling perturbative, pre- 
dicts a slight reduction of mobility in agreement with 
results present in the literature^ On the other hand, 
going from pentacene to naphthalene, the reorganization 
energy nearly increases twice suggesting a much stronger 
coupling with local modest Therefore, for systems with 
reduced number of benzene rings like naphthalene, we ex- 
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pect a larger interplay between intra- and inter-molecular 
modes within the intermediate el-ph coupling regime for 
both modes. 

In Fig. 2, we report the density of states D for 
gintra = 1-3 and Xinter — 0.09 at different temperatures. 
At T — 0, where the effect of non local modes is weak, 
the lowest band is strongly renormalized at lower ener- 
gies and becomes narrower (about half of the bare band) . 
The intrisic reduction of the bare band provides a simple 
and direct explanation of the difference in the bandwidth 
evidenced in the series of oligoacenes from naphthalene 
(effective band of the order of AOmeV) to pentacene (ef- 
fective band of the order of 80meV). Indeed, it can be 
ascribed to the decrease of the reorganization energy that 
in turn reduces the renormalization effects^ 




0.2 0.4 0.6 0.8 1.0 

T 



FIG. 3. Mobility and its different contributions as function 
of the temperature. In the inset, the mobility as function of 
temperature for different values of gintra- For comparison, 
the behavior with 5 = 2 is shown. 

At finite temperature, the shape of the spectra is 
changed due to the non local coupling. Actually, any 
band shows a new small maximum due to the coupling 
to inter-molecular modes. Moreover, the band narrow- 
ing is strongly reduced in the presence of non local cou- 
pling. In order to clarify this point, we notice that, for 
Xinter = 0, the principal band at T = 0.35 is reduced of 
about 42% of the band at T = 0. On the other hand, 
for Xinter = 0.09, the principal band at T — 0.35 is re- 
duced of only 7% of the band at T = 0. The narrowing 
of the principal band results from a subtle equilibrium 
between the two opposite tendencies. Actually, the cou- 
pling to non local modes has the main effect to induce 
scattering into the single bands of the density of states, 
preventing the narrowing induced by the coupling to lo- 
cal modes. The interplay between local and non local 
modes is able to produce a modest narrowing as function 
of the temperature even if the coupling to local modes 
is not weak. Our prediction is that this effect should be 



present not only in pentacene^ but also in naphthalene 
and anthracene. 

Finally, we analyze the mobility in the intermediate 
regime (see Fig. 3). As shown in the previous section, 
the mobility can be divided into two contributions^: the 
coherent, where the scattering of the renormalized elec- 
tron (the only effect due local el-ph coupling is here the 
reduction of the bandwidth) with non local modes is in- 
cluded, and the incoherent one, where, in addition to non 
local modes, scattering with multiple real local phonons 
is considered. The first term, relevant at low tempera- 
tures, bears a strong resemblance with the mobility of the 
system at gintra — 0, even if, as expected, it is smaller. 
The incoherent term starts at a temperature of about 
T = 0.25 and becomes predominant at higher tempera- 
tures. Actually, the coherent term only decreases with 
temperature even if the electron dynamics becomes in- 
coherent. The role of the local coupling is to promote 
an activated behavior in the incoherent regime. The co- 
operative effect is able to provide an activation energy 
A of about 0.25f, therefore less than one half of that 
for Xinter — and close to that extracted by the data 
in naphthalene^ Moreover, the upturn in the mobility 
starts at lower temperatures in comparison with the case 

Xinter = 0. 

The local coupling is able to affect but to not de- 
stroy the low temperature behavior dominated by the 
non local coupling. In the inset of Fig. 3, we report 
the low temperature mobility in logarithmic scale. Actu- 
ally, for gintra = 0, the mobility scales as 1/T 189 , while, 
with increasing gintra, the power- law becomes slightly 
less pronounced. In the case gintra — 1-3, the mobil- 
ity goes as 1/T 1 - 60 , still compatible with experiments in 
naphthalene^ 

In summary, the proposed model is able to capture 
many features of the mobility in oligoacenes. First of all, 
it has the correct order of magnitude: for T = 200AT 
(T ~ 0.2t), n ~ 1.5 hq ~ 10cm 2 /{V ■ s), while at 
T = 100K (T ~ O.li), fj, ~ 4/x ^ 28cm 2 /{V ■ s). These 
values compare favorably with mobility measurements in 
naphthalene: for T = 200K, fi ~ 4cm 2 /(V • s) in the 
hole channel and fi ~ lem 2 / (V ■ s) in the electron chan- 
nel, while at T = 100K, fx ~ 20cm 2 /(^ • s) in the hole 
channel and /i ~ 3cm 2 /(V ■ s) in the electron channel^ 
At low temperatures, the mobility scales as a power-law 
T~ s with the exponent 6 ~ 1.6 larger than unity even if 
an intermediate el-ph coupling to local modes is present. 
Starting from room temperature, the incoherent contri- 
bution of mobility begins increasing. At high temper- 
atures, the mobility shows an hopping behavior with a 
small activation energy. The theory is intended to de- 
scribe the mobility in the direction where it has the max- 
imum value. Typically, this corresponds to a direction 
into the growth plane of the OFET on the substrate^ 
The study of mobility in the other in-plane direction and 
out of plane is out of the aim of this study. 
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V. CONCLUSIONS 

The interplay between local and non local el-ph cou- 
plings strongly affects the band narrowing and the mobil- 
ity as function of the temperature. The band narrowing 
has a polaronic origin but it is strongly dependent on 
the non local el-ph coupling. In the intermediate cou- 
pling regime of the electron with both intra- and inter- 
molecular modes, the mobility scales as a power-law T~ s 



at low temperatures, and it is characterized by hopping 
behavior with a small activation energy at high temper- 
atures. We believe that the main conclusions will not be 
qualitatively modified by the inclusion of more realistic 
lattice structures and interactions since the results rely 
on general features of different types of el-ph coupling. 
The next step will be to investigate the effect of traps 
on the low temperature mobility generalizing the model 
presented in this work. 
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